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Abstract 

Using a combination of linear response theory and constrained orbital hy- 
bridization approach, we study the mechanism of magnetic exchange in- 
teraction of iron-based superconductor. We reproduce the observed highly 
anisotropic exchange interaction, and our constrain-orbital calculation un- 
ambiguously identifies that the anisotropic feature of exchange interaction is 
not sensitive to the unequal d xz /d yz orbital population. 
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The discovery of high-temperature superconductivity in iron arsenides has 
attracted intense research interests^^^. While the mediator of pairing in 
these systems remains officially unidentified, a large amount of circumstan- 
tial evidence points to magnetic spin fluctuations. Therefore, tremendous 
amount of efforts have been devoted to understand the magnetic proper- 
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However, despite vast efforts, the nature of magnetism in the iron-based 
superconductor is still a hotly debated topio^. Early theoretical studies sug- 
gest that superconducting iron arsenides have an antiferromagnetic spin- 
density-wave (SDW) instability due to Fermi-surface nesting-^. Neutron 
scattering experiment- confirms that LaFeAsO indeed exhibits the predicted 
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stripe antiferromagnetic (S-AFM) long-range ordering followed by a small 
structural distortion. However, the observed magnetic moment is much small 
than the theoretical oneA Moreover, although the general picture fits with 
a SDW model, there remain problems of matching to a purely itinerant 
scenario. In particular, the increased conductivity found in SDW state is 
not expected if a portion of the carriers become gapped 1 . Alternatively, 
a Heisenberg magnetic exchange model had been proposed to explain the 
magnetic behavior- 1 ^ 1 ^. It had been suggested that nearest-neighbor and 
next-nearest-neighbor interactions between local Fe moments are both an- 
tiferromagnetic and of comparable strength-^, which results in a magnetic 
frustration. These frustrating effects have been used to explain the structural 
phase transition and small ordered moment-. It was also suggested that the 
structural transition is actually a transition to a "nematic" ordered phase 
which will occur at a higher temperature than the SDW transition—. 

On the other hand, a short-range and highly anisotropic exchange interac- 
tion had been predicted theoretically— and confirmed by the neutron scatter- 
ing measurement subsequently 8 . To understand this unexpected anisotropy 
is a hot topio^^^. As a natural way to break the symmetry, orbital 
ordering (00) had attracted intensive research attention— liLi&ii^ and there 
is increasing experimental evidence about the orbital physics^. Band struc- 
ture calculation proposes that the degeneracy between d xz and d yz orbital 
had been lifted and there is a ferro-orbital ordering, which results in not only 
the strong anisotropic exchange but also structural transition—. However, 
the electronegativity of As is much smaller than that of O, the crystal- field 
effect upon the 3d orbitals of Fe is much weaker than in transition metal 
oxides, consequently the orbital polarization is quite small. The 00 had 
also been supported by the model calculation, but it is not clear whether the 
exchange anisotropy is related to 00 or not 17 . Therefore, a extensive study 
about the mechanism of exchange interaction is an important problem. In 
this work we address this issue using the linear response approximation— 
as well as a recently developed constrained orbital hybridization approach—. 
While our linear response approximation reproduce the known anisotropic 
exchange interaction, our constrained orbital calculation allows us to pro- 
vide theoretically a conclusive insights to various contributions to magnetic 
exchange interactions. 

We perform our electronic structure calculations based on the full-potential, 
all-electron linearized-muffn-tin-orbital (LMTO) method~. Since for this 
system local spin density approximation (LSDA) can give reasonable re- 
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Figure 1: Definitions of x, y axis, the nearest neighbor exchange interaction J% x , J\ v and 
the next nearest neighbor exchange interaction J2. 



sults^ 1 ^, we therefore adopt it as the exchange-correlation potential. With 
the electronic structure information, we estimate the exchange interaction J 
based on a magnetic force theorem— that evaluates linear response due to 
rotation of magnetic moments^. This technique has been used successfully 
for evaluating magnetic interactions in a series of compoundsJ^i'^L^. The 
main results and conclusions are found to be the same for all iron arsenides, 
we therefore focus on LaFeAsO at the following. 

The calculations are performed on the high-temperature tetragonal struc- 
ture 29 . The x and y axes are taken to be along the Fe-Fe bond direction, 
with the x axis chosen along the AFM ordered direction of S-AFM as shown 
in Fig.l. Our calculated ground state properties, including the magnetic or- 
dering configuration, density of state and band structure, are found to be in 
good agreement with previous theoretical results^. Based on the electronic 
structure information, we evaluate the interatomic exchange constants as an 
integral over the q space using (8,8,8) reciprocal lattice grid. Our numerical 
results show that despite the metallic nature, the exchange interaction is a 
short range one with the magnetic coupling further than the second nearest 
neighbor to be almost equal to zero. The short-range feature of the exchange 
interaction may be caused by the small density of state at Fermi energy. 
We reproduce the experimental observed strong anisotropic near-neighbor 
exchange interaction-. With the definition of positive J meaning the anti- 
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ferromagnetic coupling, our numerical data of J\ X1 J\ y and J 2 are 47.9, -8.0 
and 21.0 meV, respectively, which are in good agreement with the previous 
theoretical results 14 . 

Our LSDA calculation confirms that there is a small orbital polarization, 
and the difference between the occupation of d xz and d yz orbital is 0.135, 
which is very closed to the previous theoretical work (0.141) 30 . The magnetic 
moment at d xz and d yz orbital are 0.202 and 0.361 \i B respectively, which is 
also consistent with the previous calculation (0.149 and 0.338 fi B ) 30 . After 
reproducing the orbital/spin polarization, we made a calculation of J's with 
an artificial constrained external potential applied to the d xz orbital of Fe to 
adjust its energy consequently to control the orbital occupation, so that we 
can check the exact effect of unequal d xz /d yz orbital population 22 . As shown 
in Fig. 2 (a), J 2 almost does not depend on the shifting of d xz level and the 
associated orbital polarization, which is contrary to the suggestion of strong 
dependence in Ref.— . Although the value of Ji x and Ji y do depend on the 
OO, but as shown in Fig. 2(b) and Fig.2(c), even the d xy and d xz orbital has 
the same occupation (i.e., OO equal to zero), there is still strong anisotropy 
between them (J la , is almost twice larger than Ji^). Thus, we can conclude 
that the d xz and d yz orbital do have unequal population, but the anisotropic 
exchange is not related to it. 

It is well known that the strength of hybridization between two orbitals 
strongly depends on their energy difference, therefore the exchange interac- 
tion will be sensitive to the shifting of special orbital if this orbital participate 
in the exchange process. We thus perform the constrained-hybriziation ap- 
proach— to exactly analyze the possible virtual exchange mechanism directly. 
This technique has been used successfully in perovskite ruthenates and Eu- 
ropium Monochalcogenides22>2&. It turns out that a upshift of 5d orbital of 
La or a downshift of 2p orbital of O does not affect the exchange interac- 
tion. Therefore, the exchange process happens almost completely in the FeAs 
layer, and the inter-layer exchange interaction is negligible. 

In additional to d xz , we also shift other 3d orbitals of Fe. Shifting the 
3d orbitals changes the orbital occupation, however only shifting d xy orbital 
has considerable effect on J 2 . Since As anion is located above the center 
of the Fe plaquette, one can expect that the hybridization between Fe-d^ 
and As-p x ± y is strong. Thus, our numerical results clearly show that J 2 is 
mainly contributed by the As-bridged antiferromagnetic super exchange. In 
contrast to J 2 , all 3d orbitals have large effect on J la . and J\ y , which indicates 
the importance of exchange interaction due to the direct hopping between 
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Figure 2: The relation between exchange interaction and orbital polarization, (a) is J2; (b) 
is Ji x ; (c) is Ji y . The x-axis is orbital ordering (i.e. the difference between the occupation 
of d xz and d yz orbital), y-axis is the strength of exchange interaction (in meV). 
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nearest-neighbor Fe 3d electron. 

It is well known that the interatomic magnetic interaction basically is a 
band structure effect, and the spin ordering affects the covalency and details 
of the bonding topology. Therefore, it is not surprised that the exchange 
interaction depends on the magnetic configuration. For example, our ad- 
ditional calculation shows that even for NiO, which has well defined local 
moment, there is about 10% difference between the J from AFM and FM 
configuration calculation. The magnetism in iron arsenides is much more 
itinerant, moreover, there is a competition between the As-Fe superexchange 
and Fe-Fe exchange interaction. The combination of these effects results in 
the highly anisotropic nearest neighbor exchange interaction. 

To clarify the relation between the structural transition and magnetic 
property, we also perform calculation for low-temperature orthorhombic phased 
Same with the high-temperature tetragonal structure, for orthorhombic phase 
the S-AFM configuration is also lower in energy comparing with other states. 
We reproduce that the ground state is the one with the magnetic moments 
at the iron sites aligning antiparallel along the longer a axis. However, 
both the obtained magnetic moment (1.67 (j, B ) and the exchange interac- 
tion (Jia;=48.2, Jij,=-10.1, and J3=21.1 meV) are almost the same as those 
in the high-temperature phase. Moreover, we optimize lattice parameter and 
the internal atomic coordinate for both stripe antiferromagnetic ordering (S- 
AFM) and checkboard antiferromagnetic ordering (C-AFM). Our numerical 
results confirm that the structure of Fe-pnictide is almost not depend on the 
magnetic configuration. Therefore, exchange-striction effect, which had been 
used to explain the uncentrosymmetric structural distortion and the associ- 
ated multiferroics^i, cannot be used to explain the orthorhombic-tetragonal 
transition. 

In summary, based on a combination of linear response theory and con- 
strained orbital hybridization approach, we study the mechanism of magnetic 
exchange interaction of iron-based superconductor. Our results unambigu- 
ously identify that the magnetic exchange process happens in the FeAs layer, 
and the highly anisotropic feature of exchange interaction is not related to 
the orbital polarization. The magnetism is at least partially itinerant, which 
results in the anisotropic exchange interaction. While, the next nearest neigh- 
bor interaction J2 is mainly contributed by the As-bridged superexchange, 
Fe-Fe exchange interaction has considerable effect on the nearest neighbor 
exchange interaction J lx and Ji y . 

The work was supported by National Key Project for Basic Research of 



6 



China (Grant No. 2011CB922101, and 2010CB923404), NSFC under Grant 
No. 10974082, and 11174124. 

References 

[1] J. Paglione and R.L. Greene, Nature Phys. 6 (2010) 645. 

[2] M.R. Norman, Science 332 (2011) 196; F. Wang and D.H. Lee, Science 
332 (2011) 200. 

[3] LI. Mazin, Nature 464 (2010) 183; Z. P. Yin, K. Haule, G. Kotliar, 
Nature Materials 10 (2011) 932. 

[4] D.C. Johnston, Adv. Phys. 59 (2010) 803. 

[5] A.L. Wysocki, K.D. Belashchenko and V.P. Antropov, Nature Physics 
7 (2011) 485. 

[6] D. J. Singh and M.H. Du, Phys. Rev. Lett. 100 (2008) 237003; LI. 
Mazin, D.J. Singh, M.D. Johannes, and M.H. Du, Phys. Rev. Lett. 101 
(2008) 057003. 

[7] J. Dong et al, Europhys. Lett. 83 (2008) 27006. 

[8] J. Zhao et al, Nature Phys. 5 (2009) 555. 

[9] T. Yildirim, Phys. Rev. Lett. 101 (2008) 057010. 
[10] F.J. Ma, Z.Y. Lu and T. Xiang, Phys. Rev. B 79 (2008) 224517. 
[11] Q.M. Shi and E. Abrahams, Phys. Rev. Lett. 101 (2008) 076401. 
[12] C. de la Cruz et al, Nature 453 (2008) 899. 

[13] C. Fang, H. Yao, W.F. Tsai, J.P Hu and S.A. Kivelson, Phys. Rev. B 
77 (2008) 224509. 

[14] M.J. Han, Q, Yin, W.E. Pickett and SY. Savrasov, Phys. Rev. Lett. 
102 (2009) 107003. 

[15] M.D. Johannes and LI. Mazin, Phys. Rev. B 79 (2009) 220510. 

[16] C.C. Lee, W.G. Yin and W. Ku, Phys. Rev. Lett. 103 (2009) 267001. 

7 



E. Bascones, M.J. Calderon and B. Balenzuela, Phys. Rev. Lett. 104 
(2010) 227201. 

F. Kruger, S. Kumar, J. Zaanen and J. van den Brink, Phys. Rev. B 79 
(2009) 054504. 

M. Daghofer, Q.-L. Luo, R. Yu, D.X. Yao, A. Moreo and E. Dagotto, 
Phys. Rev. B 81 (2010) 180514. 

T. Shimojima et al., Phys. Rev. Lett. 104 (2010) 057002. 

X. Wan, Q. Yin and SY. Savrasov, Phys. Rev. Lett. 97 (2006) 266403. 

X. Wan, J. Zhou and J. Dong, EPL 92 (2010) 57007. 

S. Y. Savrasov, Phys. Rev. B 54 (1996) 16470. 

M. D. Johannes and I. I. Mazin, Phys. Rev. B 79 (2009) 220510. 

I. I. Mazin, M. D. Johannes, L. Boeri, K. Koepernik, D. J. Singh, Phys. 
Rev. B 78 (2008) 085104. 

A.I. Liechtenstein et al., J. Magn. Magn. Mater. 67 (1987) 65; P. Bruno, 
Phys. Rev. Lett. 90 (2003) 087205. 

X. Wan, T.A. Maier, and SY. Savrasov, Phys. Rev. B 79 (2009) 155114; 
M.J. Han, X. Wan, and S.Y. Savrasov, Phys. Rev. B 78 (2008) 060401. 

X. Wan, J. Dong and S.Y. Savrasov, Phys. Rev. B 83 (2011) 205201. 

T. Nomura, S.W. Kim, Y. Kamihara, M. Hirano, P.V. Sushko, K. Kato, 
M. Takata, A.L. Shluger and H. Hosono, Supercond. Sci. Technol. 21 
(2008) 125028. 

See the Supplemental Material of— 

S.W. Cheong and X. Mostovoy, Nature Mater. 6 (2007) 13. 



8 



